Complete Active Space multiconfiguration Dirac-Hartree-Fock calculations of 
hyperfine structure constants of the gold atom 
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The multiconfiguration Dirac-Hartree-Fock (MCDHF) model has been employed to calculate the 
expectation values for the hyperfine splittings of the 5ci^6s^ and 5d^6s^ ^Ds,/2 levels of atomic 

gold. One-, two-, and three-body electron correlation effects involving all 79 electrons have been 
included in a systematic manner. The approximation employed in this study is equivalent to a Com- 
plete Active Space (CAS) approach. Calculated electric field gradients, together with experimental 
values of the electric quadrupole hyperfine structure constants, allow us to extract a nuclear electric 
quadrupole moment Q(^^^Au)=521. 5(5.0) mb. 



PACS numbers: 31.15.am, 31.15.vj, 31.30.Gs, 21.10.Ky 



I. INTRODUCTION 



Ah initio calculations of atomic properties can now 
be performed routinely, both in the framework of the 
MCDHF theory [J, [3, y, 0,11, as well as many-body per- 
turbation (MBPT) theory Both these meth- 
ods are designed to evaluate in a systematic manner the 
electron-electron correlation effects, which constitute the 
dominant correction to all ah initio calculations based 
on the central-field approach. However, the complexity 
increases rapidly with the atomic number, and fully cor- 
related calculations, in which all electrons are explicitly 
correlated, are still possible only for very light elements 
(see e.g. [l^ [ll|, [l^l for model calculations of hyperfine 
constants of lithium- like systems) . For heavy atoms both 
theories can only be applied in a limited model (one- and 
two-body correlation effects) or only to certain atoms 
(closed-shell systems or alkali-like systems). The main 
purpose of the present paper was to carry out an ac- 
curate calculation of hyperfine structure constants of a 
heavy atom within the framework of the MCDHF the- 



ory. The calculations described in the present paper are, 
to our knowledge, the first successful evaluation of one- 
, two-, and three-body electron correlation effects for a 
heavy, open-shell, neutral atom. The multiconfiguration 
model applied in the present paper is effectively equiv- 
alent to a Complete Active Space (CAS) approach, in 
the sense that in the calculation of the hyperfine electric 
quadrupole moments all non-negligible electron correla- 
tion effects were explicitly accounted for at 1 % level 
of accuracy or better. The gold atom has been chosen. 




I, 16 1 , the nu- 

is, m 



because the hyperfine structures [Ic 
clear electric quadrupole moments 

and other properties [H, [H, HI] of gold have been the 
subject of much activity recently (the latest summary of 
nuclear quadrupole moments is given in ref. The 
second objective of the present paper is an evaluation of 
the electric quadrupole moment Q of the ^^^Au isotope. 
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II. THEORY 

The numerical-grid wavefunctions [l| were generated 
as the self-consistent solutions of the Dirac-Hartree-Fock 
equations p6l | in systematically increasing multiconfig- 
uration bases (of size NCF, which is a commonly used 
shorthand of 'Number of Configuration Functions') of 
symmetry-adapted configuration state functions $(7^ J) 



NCF 



(1) 



where ^(J) is an eigcnfunction of even parity and of 
total angular momentum J for each of the two states 
^-(5(^9652 21)3/2) and «'(5d«6s2 ^D^/^) of the isotope 
^ygAu. The sets 7^ describe multiconfiguration expan- 
sions, for which configuration mixing coefficients c^. were 
obtained through diagonalisation of the Dirac-Coulomb 
Hamiltonian 



Hoc = 



COLi ■ Pi 



V{n)] 



E 



(2) 

All calculations were done with the nucleus modelled as 
a sphere, where a two-parameter Fermi distribution (27j 
was employed to approximate the radial dependence of 
the nuclear charge density. The nuclear magnetic dipole 
moment = 0.145746(9) fi^ of ^79 Au has been used in 
calculations of magnetic dipole hyperfine constants [H, 



III. METHOD 

The numerical wave functions were obtained inde- 
pendently for the two levels of interest, 5dP6s^ ^^3/2 



and bdPGs^ 
phases: 



1)5/2. The calculations proceeded in eight 



1. Spectroscopic orbitals were obtained in the Dirac- 
Hartree-Fock approximation. These were kept 
frozen in all subsequent calculations. 

2. Virtual orbitals were generated in an approxima- 
tion (called SrD, and explained in the following 
subsection IIII A[) . in which all single and restricted 
double substitutions from SspdAspdfbspdGs spec- 
troscopic orbitals to eight layers of virtual orbitals 
were included (see the following subsection IIII Al 
for definitions of spectroscopic and virtual orbital 
sets) . 

3. Contributions from ls2sp shells were added in the 
configuration- interaction calculation, i.e., with all 
orbitals frozen. Only single substitutions con- 
tributed to the expectation values. The configu- 
rations involving ls2sp orbitals were carried over 
to the following phases. 



TABLE I: Calculated values of A and Q obtained in several 
approximations during the process of generation of virtual 
orbital set for the -D3/2 state. DHF - uncorrelated Dirac- 
Hartree-Fock value; n - largest principal quantum number 
in the orbital set; from - spectroscopic orbitals opened for 
SrD substitutions; to - virtual orbital set; NCF - number of 
configurations; see text for further details. 





experiment 




199.8425(2) 


n 


from to 


NCF 


A [MHz] Q[mh] 



DHF 


1 


218.011 


580.807 


7 


5d6s Ispdfgh 1147 


187.302 


623.275 


8 


5spd6s 2spdfgh 13729 


198.774 


652.057 


9 


4spdf...6s Sspdfgh 97526 


195.492 


547.891 


10 


3spd...6s 4spdfg3h 222129 


196.513 


528.752 


11 


3spd...6s 5spdfg3h 222494 


199.413 


523.736 


12 


3spd...6s 6spdfg3h 222851 


199.455 


514.186 


13 


3spd...6s 7spdfg3h 223212 


200.431 


515.489 


14 


3spd...6s 8spdfg3h 223573 


199.871 


515.495 



4. Unrestricted single and double substitutions (SD) 
were performed, in which one or two occupied or- 
bitals from the 5spd6s subshells were replaced by 
orbitals from the virtual set '3spdf2glh', i.e., 3 vir- 
tual orbitals of each of the 's,p,d,f ' symmetries, plus 
2 virtual orbitals of the 'g' symmetry, and 1 virtual 
orbital of the 'h' symmetry. 

5. Unrestricted triple substitutions (T) from 5spd6s 
valence and core orbitals to '2spdflg' virtual set 
were added. 

6. The final series of configuration-interaction calcu- 
lations were based on the multiconfiguration ex- 
pansions carried over and merged from all previous 
phases enumerated above. 

7. Contributions from the Breit interaction were eval- 
uated in the single-configuration approximation, in- 
cluding the full Breit operator in the self-consistent- 
field process. 

8. The values of the nuclear electric quadrupole mo- 
ment Q(^^^Au) were obtained from the relation 
B{J) = 2e(3(jJ|T(2)|Jj), where the electronic op- 
erator T(2) represents the electric field gradient at 
the nucleus. Expectation values of hyperfine con- 
stants A and of electric field gradients were cal- 
culated [3] 



5/2- 



separately for both states, and 
The experimental values of the hyperfine 



constants A and B were taken from 1311. \32 



A. Virtual orbital set 

We generated 8 layers of virtual shells (3 layers with 
'spdfgh' symmetries and 5 layers with 'spdfg' symme- 
tries). It should be noted, that the notion of a 'layer' 
is somewhat different when applied to occupied (also re- 
ferred to as spectroscopic) orbitals, as opposed to virtual 
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TABLE II: Calculated values of A and Q obtained in several 
approximations during the process of generation of virtual 
orbital set for the -D5/2 state. DHF - uncorrelated Dirac- 
Hartree-Fock value; n - largest principal quantum number 
in the orbital set; from - spectroscopic orbitals opened for 
SrD substitutions; to - virtual orbital set; NCF - number of 
configurations; see text for further details. 



40 



^ 20 - 





experiment 


80.236(3) 


n 


from to 


NCF A [MHz] Q[mb] 



DHF 


— — 1 


79.041 


612.985 


7 


5d6s Ispdfgh 11984 


69.487 


707.216 


8 


5spd6s 2spdfgh 33291 


72.278 


673.387 


9 


4spdf...6s 3spdfgh 128639 


77.761 


558.526 


10 


3spd...6s 4spdfg3h 290612 


81.020 


532.862 


11 


3spd...6s 5spdfg3h 291039 


81.045 


534.635 


12 


3spd...6s 6spdfg3h 291466 


81.248 


520.409 


13 


3spd...6s 7spdfg3h 291893 


81.214 


520.890 


14 


3spd...6s 8spdfg3h 292320 


82.136 


520.259 



(also referred to as correlation) orbitals. A core 'layer', 
i.e., a subset of occupied orbitals possessing the same 
principal quantum number (often referred to as a shell), 
constitutes a set of one-electron spin-orbitals, clustered 
in space, and having similar one-electron energy values. 
On the other hand, virtual orbitals with the same princi- 
pal quantum number are not necessarily spatially clus- 
tered because their one-electron energy values do not 
have physical meaning and may vary widely, depending 
on the correlation effects that a particular virtual orbital 
describes. Therefore a 'virtual layer' usually means a 
subset of the virtual set, generated in one step of the 
procedure, as described above. Such a 'layer' is often 
composed of orbitals with different angular symmetries. 
The notation used in the tables and text of the present 
paper reflects the above considerations, in the sense that 
occupied orbitals are listed by their principal and an- 
gular quantum numbers (i.e. 5spd means three occupied 
orbitals of s, p, and d symmetry with principal quantum 
number n = 5), while virtual orbitals are listed by an- 
gular symmetry and quantity (i.e. '5spd' would mean 
fifteen virtual orbitals — five of each of the 's', 'p', and 
'd' symmetries). To avoid confusion we distinguish oc- 
cupied orbitals from virtual ones in the present paper by 
using italics for occupied orbitals, while virtual orbitals 
are enclosed in 'quotation marks'. This distinction is not 
applied in the tables, since in the tables there are always 
headings 'from' and 'to' which clearly denote occupied 
and virtual orbitals, respectively. The notation should 
always be analysed in the proper context (see (ssj for 
further details). In the present calculations single and 
restricted double (SrD) substitutions were allowed from 
valence and core orbitals (starting from 5d6s for the first 
virtual layer). The restriction was applied to double sub- 
stitutions in such a way that only one electron was substi- 
tuted from core SspdAspdfBspd shells, the other one had 
to be substituted from valence 6s shell. Each subsequent 
layer was generated with substitutions from deeper core 
shells, down to 3s. Table[T]shows which occupied orbitals 
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5d 5p 5s 4f 4d 4p 4s 3d 3p 3s 2p 2s Is 
one-electron orbitals 

FIG. 1: Contributions from occupied orbitals to the calcu- 



lated value of Q for the 5d 6s 



■D3/2 state of Au. 



See caption 



of the Table Hn] and Sec HIlBl for further details. 



were opened at each step, as well as composition of the 
virtual orbital set when subsequent layers were generated 
for the 5d^6s^ ^^3/2 state. For instance, the line marked 
'10' in the first column describes the generation of the 
fourth virtual layer, for which the largest principal quan- 
tum number was 10; all occupied orbitals between 3s and 
6s (i.e. 3spd4:spdf5spd6s) were opened for substitutions; 
the virtual set was composed of 4 orbitals of symmetries 
's', 'p', 'd', 'f, 'g', and 3 orbitals of 'h' symmetry. 

The last four layers (those with principal quantum 
numbers 11, 12, 13, 14) were generated with a further re- 
striction, which allowed only single substitutions to these 
last layers. 

Table [II] presents the analogous data obtained in the 
process of generation of virtual orbital layers for the 
5d^6s^ state. The data from both tables arc also 

presented as crosses in figure [3] 



B. Contributions from ls2sp orbitals 

After generating the virtual orbital set, all orbitals 
were frozen and further calculations were carried out in 
the configuration-interaction (CI) approach. First, the 
effects of \s2sp orbitals were evaluated in separate CI cal- 
culations. For the 5o?^6s^ ^-03/2 state they are presented 
in the Table IIIIl together with the contributions of all 
other occupied orbitals of the gold atom. The orbitals 
that were open for single and restricted double substi- 
tutions to the full virtual set are listed in the first col- 
umn. The contributions of individual orbitals (i.e., of 
the leftmost orbital in the first column) are listed in the 
fourth and sixth column and presented in graphical form 
in Fig.[T] The individual contributions of the 2p, 2s, and 
Is orbitals to the total Q value were of the order of 0.6 %, 
0.2 %, and 0.02 %, respectively. The combined contribu- 
tion of ls2sp shells was of the order of 0.8 %, with respect 
to the total Q value. The contribution to the calculated 
value of magnetic dipole hyperfine constant A was eval- 
uated in the same manner as for Q. 
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TABLE III: Contributions from occupied orbitals to the cal- 
culated values of A and Q for the 5c?^6s'^ ^^3/2 state of Au; 
orbitals = set of orbitals open for single and restricted double 
substitutions from all shells listed in the first column to the 
full virtual set; NCF = size of the multiconfiguration expan- 
sion; AA — contribution [MHz] of the leftmost orbital from 
a given orbital set to the total A value (i.e. the individual 
contribution of the Is orbital is listed in the line Is. .6s); AQ 
= contribution [mb] of the leftmost orbital from a given set 
to the Q value. 



orbitals 


NCF 4 [MHz] 


AA 


Q[mb] 


AQ 




1 


218.011 




580.807 




5d6s 


16457 


189.406 


-28.605 


613.418 


32.611 


5pd6s 


39808 


153.103 


-36.302 


629.418 


16.000 


5spd6s 


48129 


190.849 


37.746 


625.559 


-3.859 


4f5spd6s 


89477 


187.661 


-3.188 


623.738 


-1.821 


4df5spd6s 


124673 


194.614 


6.953 


593.938 


-29.800 


4pdf5spd6s 


148188 


202.721 


8.107 


528.442 


-65.496 


4spdf5spd6s 


156525 


196.476 


-6.245 


529.646 


1.204 


3d4spdf5spd6s 


191721 


199.346 


2.870 


525.342 


-4.304 


3pd...6s 


215236 


201.106 


1.760 


514.175 


-11.167 


3spd...6s 


223573 


199.872 


-1.234 


515.495 


1.320 


2p3spd...6s 


247088 


196.564 


-3.308 


518.635 


3.140 


2sp3spd...6s 


255425 


199.576 


3.012 


519.539 


0.904 


Is2sp3spd...6s 


263762 


199.554 


-0.022 


519.634 


0.095 



A similar procedure has been carried out for the Q and 
A values of the 5d^6s^ ^D^j2 state. The results for the 
^1)5/2 state are shown in table [Wl and in figure [D The 
individual contributions of the 2p, 2s, and Is orbitals to 
the total Q value were of the order of 0.5 %, 0.2 %, and 
0.02 %, respectively. The combined contribution of ls2sp 
shells was of the order of 0.7 %, with respect to the total 
Q value. 

All these contributions have been included in the Q and 
A values obtained within the SrD approximation and the 
configuration state functions (CSFs) involved in evalua- 
tion of these contributions were carried over to all subse- 
quent calculations. 

It should be pointed out that the data in Tables IIIII 
and IIVI and in the Figures [1] and [2] were obtained with 
single and restricted double substitutions, i.e., with unre- 
stricted double and triple substitutions excluded. There- 
fore the contributions of the bpsd and Aspdf shells are 
somewhat distorted — if double and triple substitutions 
were included, the individual contributions of the 5psd 
and Aspdf shells would differ by a few percent. Only the 
3spd, 2sp and Is shells are essentially insensitive to dou- 
ble and triple substitutions (see Sec. lIII Cl below). There- 
fore their contributions are approximately correct. 



C. Double, triple, and quadruple substitutions 

The decomposition of the electron correlation correc- 
tion to the hyperfine structure into one-, two-, three-, 
and four-body effects can be understood from the follow- 
ing (simplified) analysis. The structure of the 5d^6s^ 
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5d 5p 5s 4f 4d 4p 4s 3d 3p 3s 2p 2s Is 
one-electron orbitals 

FIG. 2: Contributions from occupied orbitals to the calcu- 
lated value of Q for the 5d''6s^ ^^5/2 state of Au. See caption 
of Table HV] and Sec HITBl for further details. 



TABLE IV: Contributions from occupied orbitals to the cal- 
culated values of A and Q for the 5d^6s^ ^^^5/2 state of Au; 
orbitals = set of orbitals open for single and restricted double 
substitutions from all shells listed in the first column to the 
full virtual set; NCF = size of the multiconfiguration expan- 
sion; AA = contribution [MHz] of the leftmost orbital from 
a given orbital set to the total A value (i.e. the individual 
contribution of the Is orbital is listed in the line Is. .6s); AQ 
= contribution [mb] of the leftmost orbital from a given set 
to the Q value. 



orbitals 


NCF 


A [MHz] 


AA 


Q[mb] 


AQ 




1 


79.041 




612.985 




5d6s 


21501 


106.724 


27.683 


651.547 


38.562 


5pd6s 


51800 


109.554 


2.830 


643.451 


-8.096 


5spd6s 


62536 


71.472 


-38.082 


638.694 


-4.757 


4f5spd6s 


117626 


70.636 


-0.836 


636.280 


-2.414 


4df5spd6s 


163739 


73.490 


2.854 


604.280 


-32.000 


4pdf5spd6s 


194221 


74.597 


1.107 


532.632 


-71.648 


4spdf5spd6s 


204973 


79.767 


5.170 


534.008 


1.376 


3d4spdf5spd6s 


251086 


80.877 


1.110 


530.308 


-3.700 


3pd...6s 


281568 


80.580 


-0.297 


518.734 


-11.574 


3spd...6s 


292320 


82.136 


1.556 


520.259 


1.525 


2p3spd...6s 


322802 


81.700 


-0.436 


522.677 


2.418 


2sp3spd...6s 


333554 


78.995 


-2.705 


523.757 


1.080 


Is2sp3spd...6s 


344306 


79.025 


0.030 


523.880 


0.123 



states of gold is determined to a large extent by the inter- 
action of the valence 6s^ shell with a highly polarisable 
5d^ shell. The direct and indirect effects of relativity 
bring the outer d shell much closer, radially and ener- 
getically, to the valence s orbital than in homologous 
silver and copper atoms [13, [l^l ■ This in turn increases 
the polarisation of the 5d^ shell by the valence electrons. 
Therefore, the core- valence interaction (the leading elec- 
tron correlation correction) leads to the contraction of 
the 6s orbital, which overestimates the hyperfine struc- 
ture. The unrestricted double substitutions affect the 
hyperfine structure in two ways: directly through the 
configuration state functions (CSFs) themselves but also 
indirectly through the change of the expansion coeffi- 
cients of the important configurations obtained by sin- 
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TABLE V: Values of Q and A for the -D3/2 state, calculated 
in configuration-interaction approach, with single and unre- 
stricted double substitutions, in several different multiconfig- 
uration expansions; from - spectroscopic orbitals opened for 
substitutions; to - virtual orbital set; NCF - number of con- 
figurations; see text for further details. 



experiment 




199.8425(2) 




from 


to 


NCF 


.4 [MHz] 


Q[mb] 


5spd6s 


Ispdfgh 


259135 


205.426 


521.191 


4spdf5spd6s 


Ispdfgh 


358019 


205.968 


521.503 


5spd6s 


2spdf 


279559 


210.523 


509.839 


5spd6s 


2spdfg 


320545 


211.512 


511.461 


5spd6s 


2spdfgh 


366257 


211.480 


512.286 


5spd6s 


3spdf2glh 459594 


213.088 


510.451 


5spd6s 


3spdf2gh 


465794 


213.075 


510.402 


5spd6s 


3spdfg2h 


506987 


213.146 


510.268 


5spd6s 


4spdf2gh 


687301 


213.200 


510.478 



TABLE VIL Values of Q and A for the -D3/2 state, calculated 
in configuration-interaction approach, with single and unre- 
stricted double and triple substitutions, in several different 
multiconfiguration expansions; from - spectroscopic orbitals 
opened for substitutions; to - virtual orbital set; NCF - num- 
ber of configurations; see text for further details. 



experiment 


199.8425(2) 




from 


to 


NCF 


A [MHz] 


0[mb] 


5spd6s 


Ispd 


265183 


198.955 


520.346 


5spd6s 


Ispdf 


386326 


194.391 


533.464 


5spd6s 


Ispdfg 


641227 


193.744 


536.620 


5spd6s 


Ispdfgh 1012615 


193.246 


537.786 


5spd6s 


2spdlf 


943544 


198.752 


522.361 


5spd6s 


2spdf 


1543051 


199.973 


520.536 


5spd6s 


3spd2f 


1200261 


198.207 


520.267 


5spd6s 


3psdf 


1309130 


198.254 


520.096 



TABLE VI: Values of Q and A for the -D5/2 state, calculated 
in configuration-interaction approach, with single and unre- 
stricted double substitutions, in several different multiconfig- 
uration expansions; from - spectroscopic orbitals opened for 
substitutions; to - virtual orbital set; NCF - number of con- 
figurations; see text for further details. 



experiment 




80.236(3) 




from to 


NCF 


A [MHz] 


Q[mb] 


5spd6s Ispdfgh 


339306 


74.258 


507.823 


4spdf5spd6s Ispdfgh 


467381 


72.048 


509.321 


5spd6s 2spdfgh 


480824 


73.468 


512.278 


5spd6s 3spdf2gh 


607421 


73.494 


512.559 


5spd6s 4spdf3gh 


898368 


73.294 


514.621 


5spd6s 5spdf4g3h 1228675 


73.212 


514.269 



gle substitutions. Three-particle effects in turn affect 
the expansion coefficients of the configurations obtained 
from double substitutions. In a simple picture we can 
describe the wave function in terms of pair-correlation 
functions and the three-particle effects then account for 
polarisation of pair-correlation functions, leading to an 
increase of the hyperfine structure [11] . Four-particle ef- 
fects affect mostly the expansion coefficients of the config- 
urations obtained from double substitutions. Therefore 
their influence on the hyperfine structure is indirect and 
second-order to that of the double substitutions. They 
are usually small and can often be neglected [s^; they 
are discussed in Sec. IIII CI 

Tables |V] and IVII show the results of configuration- 
interaction calculations, where various combinations of 
occupied and virtual sets were tested with single and 
unrestricted double substitutions. The data from both 
tables are also presented as empty circles in Fig. [3] 

The second line in Tables IVl and IVII represents a calcu- 
lation in which substitutions from the Aspdf shells were 
allowed to one layer of virtual orbitals. When compared 
with the first line, it yields the effect of Aspdf shells on the 
calculated values of Q and A. In order to limit the size of 
the configuration expansions, the CSFs representing the 
above substitutions were not carried over to the foUow- 



TABLE VIII: Values of Q and A for the -D5/2 state, calculated 
in configuration-interaction approach, with single and unre- 
stricted double and triple substitutions, in several different 
multiconfiguration expansions; from - spectroscopic orbitals 
opened for substitutions; to - virtual orbital set; NCF - num- 
ber of configurations; see text for further details. 



experiment 


80.236 (3) 




from 


to 


NCF 


^[MHz] 


Q[mb] 


5spd6s 


Ispd 


341440 


81.6955 


514.929 


5spd6s 


Ispdf 


456506 


82.1357 


520.259 


4f5spd6s 


Ispdf 


1403860 


79.9343 


518.380 


5spd6s 


Ispdfg 


842883 


80.2371 


519.291 


5spd6s 


2spdf 


1326851 


83.0623 


521.862 



ing, higher order calculations. Instead, the corrections 
were included additively, as described in Sec. IIV Al At 
the same time the evaluation of these corrections may be 
treated as a crude estimate of error arising from omitted 
double substitutions from occupied shells (see Sec. IIVBI 
for details). 

An inspection of the last column of Tabic \V\ indicates, 
that three layers of virtual orbitals were necessary to 
reach convergence of the Q and A values in the single 
and unrestricted double substitutions (SD) approxima- 
tion for the 5d^6s^ ^^3/2 state. Four layers were neces- 



TABLE IX: Values of Q and A for the -D3/2 state, calcu- 
lated in configuration-interaction approach, with single and 
unrestricted double, triple, and quadruple substitutions, in 
several different multiconfiguration expansions; type - sub- 
stitution multiplicity; from ~ spectroscopic orbitals opened 
for substitutions; to - virtual orbital set; NCF - number of 
configurations; see text for further details. 



experiment 



199.8425(2) 



type from 



to 



NCF A [MHz] Q[mb] 



SDT 5spd6s Ispd 386326 194.391 533.464 

SDTQ 5spd6s Ispd 569497 194.301 533.653 

SDT 5spd6s Ispdf 386326 194.391 533.464 

SDTQ 5spd6s Ispdf 967871 195.376 531.685 

SDTQ 5spd6s Ispdf 1089014 194.686 531.846 
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sary in case of the Sd^Gs^ ^^5/2 state (see Tabic IVH) . 

Tables IVlIl and IVIIII show the results of configuration- 
interaction calculations, in which various combinations of 
occupied and virtual sets were tested with unrestricted 
double and triple substitutions. The data from both ta- 
bles are also presented as full circles in Fig. [3l Two layers 
of virtual orbitals were necessary to reach convergence of 
the Q value in the single and double and triple substi- 
tutions (SDT) approximation for both 5(P6s^ ^-03/2 and 
5(i^6s^ states. In case of the A values, convergence 

required three, rather than two, layers. 

Table IIXI shows the effect of quadruple substitutions. 
The first line represents an approximation in which sin- 
gle, double, and triple substitutions from 5spd6s orbitals 
to a truncated virtual layer composed of 's', 'p', and 'd' 
symmetries were included. The third line represents a 
similar approximation in which the (still truncated) vir- 
tual layer was composed of 's', 'p', 'd', and 'f symmetries. 
The second, fourth, and fifth lines represent correspond- 
ing 'quadruple approximation' in which single, double, 
triple, and quadruple substitutions were allowed. The 
comparison had to be made on a reasonably small or- 
bital set in order to be able to converge the calculation 
involving quadruple substitutions. The numbers of CSFs 
in the last two lines arc different, because certain restric- 
tions were applied in the calculation represented by the 
fourth line (see the comments near the end of Sec. IIII Dl 
for details) . The results presented in the table llXl indicate 
that the correction involving quadruple substitutions is 
unlikely to exceed 1 %. The CSFs representing quadru- 
ple substitutions were not carried over to the following 
calculations and the 'quadruple' correction was included 
additively, as described in Sec. IIV Al 



D. 4-D configuration-interaction calculations 

A full, converged Complete Active Space (CAS) cal- 
culation for the gold atom is still unattainable due to 
software and hardware limitations. Based on our cur- 
rent calculations wc estimate that the CAS approach 
would require configuration expansions in four 'dimen- 
sions': (1) single, double, triple, and perhaps quadruple 
substitutions; (2) from all core shells (or at least from 
ispdAspdfbspdQs)] (3) to eight or more virtual orbital 
layers; (4) of s, p, d, /, g, h, and perhaps higher sym- 
metries. One can imagine a 'space' spanned by the four 
'dimensions' defined above, i.e. ^substitution multiplic- 
ity', 'number of opened core subshcUs', 'number of vir- 
tual layers', and 'maximal symmetry of virtual layers' 
dimension. In fact, this space should rather be called 
a 'matrix', since all four dimensions are discrete. Let's 
call this four-dimensional matrix a 'CAS matrix'. Each 
element of the matrix is represented by a multiconfig- 
uration expansion obtained by substituting a particular 
number of electrons ('substitution' dimension) from spe- 
cific core orbitals ('core' dimension) to a set of virtual or- 
bitals ('virtual' dimension) of specific symmetries ('sym- 



TABLE X: The final configuration-interaction calculations of 
Q and A for the 5d^6s^ ^£'3/2 state of Au; type — description 
of the multiconfiguration expansions, see text for details; NCF 



size of the multiconfi, 


guration 


expansion. 




experiment 




199.8425(2) 




type 


NCF 


yl[MHz] 


Q[mb] 


SD:3hgg+ SDT:2fd 


1182329 


206.343 


517.201 


SD:3hgf H- SDT:2fd 


1144532 


206.221 


517.342 


SD:3hgf -\- SDT:2gd 1711382 


205.104 


519.106 


SD:3hgf + SDT:2gf 


1847380 


204.489 


519.829 



TABLE XI: The final configuration-interaction calculations of 
Q and A for the 5d^6s^ ^^5/2 state of Au; type — description 
of the multiconfiguration expansions, see text for details; NCF 



size of the multiconfiguration 


expansion. 




experiment 


80.236 (3) 




type NCF 


yl[MHz] 


Q[mb] 


SD:3hgf + SDT:2fd 1441120 


78.2451 


520.073 


SD:3hgf + SDT:2gd 1527668 


79.9182 


522.066 



metry' dimension). A full CAS calculation would require 
several orders of magnitude larger configuration expan- 
sions than are possible even with the largest computer 
resources available today. 

However, a computational strategy can be designed in 
which a considerably smaller multiconfiguration expan- 
sion yields a wave function only marginally inferior to a 
full CAS wave function, in the sense that all important 
electron correlation effects are included and the calcu- 
lated values of A and Q are close to those that would 
result from a full, converged, CAS calculation. The strat- 
egy is based on the observation that one does not have to 
simultaneously push the configuration expansions to the 
limits of all the above mentioned 'dimensions'. Specifi- 
cally, the dependence of atomic properties on the 'substi- 
tution' dimension is critical. To illustrate this approach, 
let's consider separately the contributions of single, dou- 
ble, and triple substitutions to the calculated values of A 
and Q of gold. To obtain converged results within a sin- 
gle substitutions model, one has to include substitutions 
from all occupied shells {ls2spSspd4:spdf5spd6s) to eight 
or more virtual layers. This is illustrated in tables |T1 [Til 
mil and IIVI where eight virtual orbital layers were nec- 
essary to converge the series of self-consistent-field cal- 
culations. However, to obtain a converged result within 
a single and double substitutions model (SD) one has to 
include double substitutions from Aspdf^spdGs occupied 
orbitals, not to eight, but to three or at most four virtual 
layers (see Tables IVl and IVip . In the single, double, and 
triple substitutions model (SDT) it is enough to consider 
triple substitutions from 5spd6s occupied orbitals to two 
or at most three virtual layers (see Tables IVlIl and IVIIip . 
In the 'space' (or rather in the 'matrix') of the four 'di- 
mensions' defined above, the 'core' and 'virtual' dimen- 
sion sizes strongly depend on the 'substitution' dimension 
(in fact, all four dimensions are interdependent). 

Therefore, one can construct an approximation, in 
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which all important electron correlation effects are in- 
cluded and the calculated values of A and Q are close to 
those that would result from a full, converged CAS calcu- 
lation. In order to find a suitable approximation, we have 
performed a set of test calculations for several elements of 
the above mentioned 'matrix'. For each 'dimension', the 
calculations were saturated to the point where the rela- 
tive change of the expectation values (i.e., both A and Q) 
did not exceed a small fraction of a percent (usually two 
or three tenths of a percent) . Specifically, for each 'sub- 
stitution' dimension (i.e., for single, double, and triple 
substitutions) we thoroughly tested the dependence of 
observables on 'symmetry,' 'virtual', and 'core' spaces. 
When a saturated set of configuration state functions 
(CSFs) is obtained for a particular 'substitution' dimen- 
sion, all these CSFs are carried over to the next stcp(s). 
The merged, 'final' multiconfiguration expansion repre- 
sents an approximation, which is effectively equivalent to 
a CAS expansion, and the corresponding wavefunction is 
of similar quality as a CAS wavefunction, at least from 
the point of view of the calculated values of A and Q. 

In practice there is not one single final 'CAS' expan- 
sion, but a series of such final expansions in which various 
sets of 'S', 'SD', and 'SDT' multiconfiguration expansions 
(i.e. various sets with single, double, and triple substitu- 
tions) are merged together. Table |X] shows the results 
obtained from a series of such final 'CAS' calculations 
for the ^-D3/2 state, and Table IXj shows the same for the 
^1)5/2 state. The data from both tables are also included 
in Fig. [3l The 'CAS' expansions are composed as fol- 
lows. All virtual orbitals and all CSFs generated in the 
SrD approximation, as described in section fill Ai as well 
as those described in section IIII B[ were included. The 
remaining CSF expansions were generated with substitu- 
tions from 5spd6s orbitals to virtual sets described in the 
first column of Tables |X] and IXIl where symbols before 
the colon represent substitution multiplicity, i.e., SD — 
single and double substitutions; SDT — single and dou- 
ble and triple substitutions; the symbols after the colon 
represent virtual orbital layers, i.e., 3hgg — three lay- 
ers (first layer with 'spdfgh' symmetries and two layers 
with 'spdfg' symmetries); 2fd — two layers (first layer 
with 'spdf symmetries and second layer with 'spd' sym- 
metries); 3hgf — three layers (first layer with 'spdfgh' 
symmetries, second layer with 'spdfg' symmetries; third 
layer with 'spdf symmetries), etc. 

In the largest calculations, when single, double, and 
triple substitutions to two or three layers were included, 
we had to further limit the overall number of CSFs, due 
to software and hardware limitations. In those cases, the 
occupation number of the least important virtual orbital 
was restricted to single or double, thus excluding those 
CSFs in which this particular virtual orbital was occupied 
by three electrons. The difference that such a restriction 
brings about can always be evaluated on a smaller set of 
CSFs before a full calculation is performed. Therefore we 
always had control on the effects of the above mentioned 
restrictions on the calculated values of A and Q. 



IV. RESULTS 

More extensive calculations turned out to be beyond 
the 100 node limit for this project on the Linux clus- 
ter at the National Institute of Standards and Tech- 
nology (NIST), USA. Therefore the calculations of the 
magnetic dipole constants A did not yield converged re- 
sults. As might be expected, the effects of double and 
triple substitutions are relatively larger for A than for 
Q, therefore the calculations of the Q values were essen- 
tially converged; they yield: Q(^£'3/2) = 519.829 mb, and 
Q(^^5/2) = 522.066 mb, respectively. 



A. Corrections 

As mentioned in section IIII C[ the contributions aris- 
ing from unrestricted double substitutions from Aspdf 
orbitals were evaluated separately and included addi- 
tively in the final Q values. They yield -1-0.312 mb and 
4-1.498 mb for the two states ^I?3/2 and ^1)5/2, respec- 
tively. The effects of the quadruple substitutions were 
also evaluated separately, in a very limited fashion, and 
only for the ^I?3/2 state. As explained in section UlI C[ 
the correction arising from the quadruple substitutions 
for the ^-03/2 state lowers the Q value by 1.779 mb. The 
dependence of the Q values on double and triple sub- 
stitutions indicates that the quadruple correction might 
be smaller for the value than for the Q(^D^/2) 

value, but we were unable to evaluate the former. There- 
fore we assumed identical, —1.8 mb corrections for both 
states. The corrections arising from the Breit interaction 
were calculated at the Dirac-Hartrec-Fock level with full 
relaxation i.e. with a frequency-dependent Breit term 

Bij = ^ — ^ — i[cos(w.yry) - 1] 

+c^(a..V.)(a,.V,) '^°^^"'^^'^/^)-^ (3) 

included in the self-consistent-ficld functional, using the 
mcdfgme code [H, 0, [s^. In the formula above, = 
\ri — rj\ is the inter-electronic distance, ujij is the en- 
ergy of the photon exchanged between two electrons, 
cx-i are Dirac matrices, and c = 1/a is the speed of 
light. The Breit corrections are highly state-dependent 
(see also [l^, where the Gaunt part was evaluated) 
and yield 2.3 mb and 0.6 mb for the two states, ^1)3/2 
and ^-D5/2, respectively. The quantum electrodynam- 
ics (QED) corrections to the Q values are expected 
to be very small. We evaluated the VP (vacuum po- 
larisation) correction with the mcdfgme code, follow- 
ing Ref. [i^l, and obtained a value of the order of 
0.01 %. When all above mentioned corrections are in- 
cluded, the Q values become: Q(^D^/2) = 520.641 mb 
and Q{^Dt^/2) = 522.364 mb. The average of the above 
two results yields Q(^^^Au) = 521.5 mb. 
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FIG. 3: (Color online) Nuclear quadrupole moment Q(^^'^Au) [mb] obtained from the calculated electric field gradients, and 
liyperfine magnetic dipole constants A [MHz] of the states 5d^6s^ and Sd^Bs^ '^^5/2, as functions of the size of the multi- 

configuration expansions; line with no symbols (red online) - 'SrD' approximation; circles (blue online) - 'SD' approximation; tri- 
angles (black online) - 'SDT' and final 'CAS' approximations (see text for details). Horizontal straight lines in figures (c) and (d) 
represent the experimental values of hyperfine constants Ai^D^ri) = 199.8425(2) MHz [H and ^(^Ds/a) = 80.236(3) MHz [sl], 
respectively. The small corrections described in subsection IIV Al are not included in the figures. 



B. Error estimate 

A rigorous, systematic treatment of the error bar of the 
calculated electric quadrupole moment Q would require 
evaluation of the effects of: all omitted virtual orbitals, 
all CSFs which were not included in the configuration 
expansions, as well as all physical effects that were not 
included or were treated approximately. However, we 
were only able to obtain very crude estimates of certain 
sources of systematic errors. We believe that none ex- 
ceeded 1 %, but the calculations presented in this paper 
were far too extensive to permit a rigorous treatment of 
the error. Therefore we have to resort to a less rigorous 
method. 

One of the frequently used methods of evaluation of 
the accuracy of calculated electric quadrupole moments 



Q is based on the simultaneous calculations of magnetic 
dipole hyperfine constants A, and on subsequent compar- 
ison of calculated A values with their experimental coun- 
terparts. As mentioned above, the calculations of the 
magnetic dipole constants A have not converged. How- 
ever, the amplitudes of the final oscillations of the two 
curves representing the values of A for the two states of 
interest are comparable to the uncertainty of A arising 
from the accuracy of the nuclear magnetic dipole moment 
value 

There are currently two different /i values in the litera- 
ture m,!!^, = 0.145746(9) and = 0.148158(8), which 
differ by about 2 %. Taken at face value, our results seem 
to favor the smaller value, = 0.145746(9), which, as 
mentioned in Sec. |TT1 has been used in the present calcu- 
lations. However, the overall accuracy of our calculations 
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size of MCDHF space 



FIG. 4: (Color online) Nuclear quadrupole moment Q [mb] of 
the ^^^Au isotope obtained from the calculated electric field 
gradients, as a function of the size of multiconfiguration ex- 
pansions for the states 5d^6s^ ^0^/2 (triangles - red online) 
and bd?Qs^ ^£^5/2 (circles - blue online), compared with other 
theoretical and experimental results. The small corrections 
described in subsection IIV Al are not included in the figure. 
The values which represent multiconfiguration expansions of 
sizes smaller than 20000 are outside the figure, except the un- 
correlated Dirac-Hartree-Fock values, represented by the sin- 
gle triangle (red online) for ^1)3/2 and the single circle (blue 
online) for The six values with error bars are from ref- 

erences p^ . [T9I . [20I . |3^ . I37I . |38| ; the four values without error 
bars, represented by pluses, are from references [itI, [21I. [sH : 
all data are arranged in reverse chronological order, with the 
most recent results to the left. 



(in particular, the evaluation of higher-order terms) does 
not permit us to draw a definitive conclusion. Therefore, 
the difference between the two values of fji should rather 
be treated as a source of systematic error in the determi- 
nation of A. Therefore, we did not push the calculations 
of magnetic dipole constants A further beyond their cur- 
rent level of convergence and, consequently, the calcula- 
tions of A values could not be used as reliable sources of 
error estimate for nuclear moments. 

Another method to estimate the accuracy of Q is to 
consider the differences between the final values obtained 
from different states. However, in the present paper we 
were able to converge the calculations for only two atomic 
levels. The difference between the results obtained for 
these two levels turned out to be quite small, which ren- 
dered this method useless in this particular case. 

Considering the computational methodology employed 
in this paper, it is obvious that the final value depends 
on the choice of the multiconfiguration expansions repre- 
senting the last few points on the curves in Fig. [H while 
the accuracy of the final value is connected with conver- 
gence of these curves. Therefore, we based the estimate 
of the error bar on the oscillations of the tail of the two 
curves in Fig.[51 The largest difference taken from the last 
few points on the curves representing ^Dj,/2 and ^I?5/2 



states amounts to 3 mb and 4 mb, respectively. As an 
additional source of uncertainty we assumed the additive 
corrections described in subsection llV Al since all of them 
were evaluated in a rather crude approximation. For in- 
stance, the contribution of the Breit interaction was cal- 
culated at the Dirac-Hartree-Fock level, without regard 
for electron correlation effects. When all above sources 
of uncertainly are taken into account the total error bar 
amounts to 5 mb, which yields our final calculated value 
of quadrupole moment Q(^'^^Au) = 521.5 ± 5.0 mb. 



V. COMPARISONS 

The results of our calculation are compared with pre- 
vious evaluations in Table IXlIl and in Fig. [4] It is worth 
noting that our result is in agreement with three most 
recent theoretical values, obtained with three different 
methods, but all these recent results (including ours) are 
considerably smaller than other, earlier values. 

Yakobi et al [l^ performed calculations for the 
bdPQs^ and 5d^6s^ ^^5/2 states of atomic 

gold within the four-component Dirac-Coulomb frame- 
work [ill, ll^l. They correlated 51 out of the 79 electrons 
in the large basis sets (up to 26s22pl8dl2f8g5h uncon- 
tracted Gaussian functions) with the relativistic Fock- 
space coupled cluster method, including single and dou- 
ble excitations (CCSD). The contribution of the Gaunt 
term, the main part of the Breit interaction, was also 
evaluated. 

Belpassi et al [2^ performed molecular relativistic 
Dirac-Coulomb-Gaunt Hartree-Fock calculations [i^ for 
a series of molecules: AuF, XeAuF, KrAuF, ArAuF, 
(OC)AuF, and AuH. The electronic correlation contri- 
butions were included at CCSD(T) and CCSD-T lev- 
els. The value of the nuclear quadrupole moment Q 
was obtained from the determinations of the electric field 
gradient at the gold nucleus for the above mentioned 
molecules, combined with experimental values of the nu- 
clear quadrupole coupling constants. 

Thierfelder et al plj performed four-component rel- 
ativistic density-functional (DFT) calculations for di- 
atomic compounds CuX and AuX (X = H, F, CI, Br, 
and I) with and without CO attached, i.e., OC-CuX and 
OC-AuX (X = F, CI, Br^ and I). They employed a newly 
developed functional [4J|, whose role is to correctly de- 
scribe the long-range part of exchange interactions [i^ , 
and obtained the averaged result Q = 526 mb. This value 
is within the error bounds of our value. 

Our result, in turn, falls within the error bounds pub- 
lished by Belpassi et al (Q ~ 510(15) mb), as well as those 
by Yakobi et al (Q = 521(7) mb). The agreement with 
Yakobi et al may be somewhat accidental because partic- 
ular contributions show larger differences. The two out- 
standing differences arise from triple substitutions and 
from deep core orbitals. Yakobi ct al evaluated the effect 
of the triple substitutions by performing single-reference 
CCSD(T) calculation for the level, and obtained 
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TABLE XII: Comparison of the present QC^'^Au) value 
(in mb = lO"'^^ m^) with other recent values and with previ- 
ous (muonic) standard value. 



Reference 



Source 



Q(^'^^Au) 



An atom, D 



3/2, 



D5/2 521.5±5.0 



An atom, ^1)3/2, ^1)5/2 521±7 



This work 
Yakobi et al [19] 

Belpassi et al [20] AuF, LAuF molecules 510±15 

Thierfelder et al [21J AuX, LAuX molecules 526 

Powers et al [38] Muonic 547±16 



a 0.3 % shift. The effect of triple substitutions is indeed 
smaller for the ^1)5/2 level, but for the ^I?3/2 level our 
calculations indicate a shift of the order of 2 %. However, 
this discrepancy may be attributed to the methodological 
differences in the two papers. The definition of triple sub- 
stitutions in the configuration-interaction method used 
here differs substantially from that in the CCSD(T) ap- 
proach, due to the exponential nature of the coupled- 
cluster operator. The coupled-cluster approximation in- 
cludes a subset of the CI triple substitutions (the "un- 
linked" diagrams), as well as that of higher order substi- 
tutions, aheady at the CCSD level. The CCSD(T) yields 
only the "linked" part as the effect of the triple substi- 
tutions. Therefore, the contribution of the CI triple sub- 
stitutions may indeed be expected to be larger than that 
of the CC triple substitutions. 

Another difference arises from contributions of deep 
core orbitals. The effects of 3spd, 2sp, and Is orbitals 
were neglected by Yakobi et al, while in our calculations 
they were all included. Their combined effect was to 
lower the Q value by about 2 %. 



Conclusions 

The multiconfiguration Dirac-Hartree-Fock (MCDHF) 
model has been employed to calculate the expectation 
values responsible for the hyperfine splittings of the 
5d^6s^ ^0^/2 and bd^Gs'^ levels of atomic gold. To 

our knowledge, this is the first calculation in which one-. 



two-, and three-body electron correlation effects were in- 
cluded and saturated for electric quadrupole hyperfine 
values of a heavy, open-shell, neutral atom. The correla- 
tion effects involving all 79 electrons were accounted for 
with a procedure that is equivalent to a full Complete 
Active Space calculation. All electron correlation effects 
were explicitly accounted for at a 1 % level of precision or 
better. Calculated electric field gradients, together with 
experimental values of the electric quadrupole hyperfine 
structure constants B, allow us to extract a nuclear elec- 
tric quadrupole moment Q = 521.5(5.0) mb of ^^^Au. If 
taken at face value, the summary in table IXIII sug gest s 
that our Q value, together with that of Yakobi et al jl9j . 
could be the new standard value. 
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